Geometrical pumping in spin coupled double quantum dots 
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We analytically investigate the non-equilibrium pumping for a double quantum dots system on 
the basis of the quantum master equation (QME), where the double quantum dots are connected 
to two external leads by the spin coupling. Each of leads has two tunable parameters, temperature 
and chemical potential. Using QME formalism, we obtain an analytical expression of the Berry-like 
phase for eigenstate of the QME in the parameter space. We show that the Berry-like curvature is 
non zero in the whole region in the parameter space. We also show that the Berry-like curvature 
vanishes in the case of single quantum dot and the case of isolated two dots. 
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I. INTRODUCTION 

Quantum transport phenomenon has been investigated 
from the early stage of the construction of the quantum 
theory. It still has attracted theoretical and experimental 
investigations to reveal the quantum many-body proper- 
ties under various non-equilibrium circumstances. 

The recent development of the nanotechnology allows 
us to make an artificial atom with nanometer scale, so 
called quantum dots, in semiconductor devices. There 
are a number of tunable parameters in the quantum dots, 
e.g., energy level in a quantum dot, bias voltage, and 
tunnel barrier between the quantum dots and the leads. 
Those tunablities enable us to compare the theory with 
the experiment for a variety of physical effects. 

One of the most cerebrated development for the quan- 
tum transport is the observation of the quantum pump in 
the quantum dot system. The original idea of the quan- 
tum pumping has been proposed by Thoulessi When 
some parameters of the Hamiltonian are adiabatically 
varied, the eigenstates obtain the Berry phase, which 
depends on the trajectory in the parameter space and 
behaves as a vector fields If the system offers non-zero 
curvature of the "vector field" in the parameter space, 
the cyclic modulation of the parameter yields currents 
originated from the topology of the system. The origi- 
nal idea of the quantum pumping for closed system has 
been extended to an open system^— and represented by 
the geometrical expression^ Since then, the various ef- 
fects on the quantum pumping have been investigated, 
e.g. Rabi oscillation between states of a double quan- 
tum dot^ interaction of two electrons in a triple- well 
structure^ Kondo effect in the Toulouse limit^ and 
so on. Experimentally, the quantum pumping has been 
realized by the transport experiment in the mesoscopic 
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In those experiments, the quantized dc cur- 



rent has been obtained in the absence of any external bias 
voltage. 

The topological phase also appears in the Master equa- 
tion. The parameter of the Liouvillian is adiabatically 



varied, the eigenstates obtain the topological phase sim- 
ilar to the Berry phase. This topological phase, so called 
Berry-Sinitsyn-Nemenman (BSN) phase, has originally 
offered in the context of the classical master equation for 
stochastic systems i 13 ' — In the systems, the cumulant 
generating function of the pumped current is expressed 
by a Berry-like phase on the eigenstate of the classical 
master equation. 

Recently, the quantum adiabatic pumping on the 
basis of quantum master equation (QME) has been 
developed^ In Ref. HH, they have analyzed the QME for 
adiabatic modulation of the reservoir parameters, such 
as temperatures and chemical potentials in leads. They 
have demonstrated that the BSN phase exists in general 
in the parameter space, and thus, the pumping current 
can exist in general situations under the adiabatic modu- 
lations of reservoir parameters. They have applied their 
method to a double quantum dots system with inter-dot 
repulsion. They obtained the BSN curvature for various 
interaction strength and have shown that the curvature 
vanishes in no-interaction limit. It is remarkable that 
the quantum pumping is generated by the modulation of 
the parameters not for systems but for thermal baths. 
However, their calculation is numerical and explicit pa- 
rameter dependence is not clear. Moreover, they have 
not taken into account the spin degrees of freedom in 
the quantum dots nor the precise intra-dot interaction 
through localized spins for double dots. The BSN cur- 
vature in their model has peaks around energy levels of 
the dots, and that implies the charge fluctuation assists 
the quantum pumping. In this paper, we analyze a more 
realistic system for double quantum dots. We calculate 
the BSN curvature analytically. We have shown that the 
spin degrees of freedom, or spin fluctuation also play an 
important role in the quantum pumping. 

The present paper is organized as follows. In Sec. II, we 
introduce the model for the double quantum dot coupled 
to external leads with spin coupling. In Sec. Ill, we make 
a brief review of calculation for the pumped current in 
the basis of quantum master equation. In Sec. IV, we 



apply the method presented in the Sec. Ill to the double 
quantum dot system. We derive the analytical expression 
for the BSN curvature in Sec. IV A in leading order of 
the perturbation. In Sec. IV B, we calculate the pumped 
current. Finally, in Sec. V we summarize and discuss our 
results. 
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II. MODEL 

In this section we present our model for a system of 
double quantum dots connected to the source and drain 
leads. In order to clarify the effect of the spin fluctuation 
on the quantum pumping, we focus on the spin local- 
ized region in which a number of the electrons in each 
quantum dot is fixed to one. In the region, the system is 
expressed by the Kondo model where the localized spin 
in the quantum dots couples to the spins in the Fermi 
see of the leads. We assume that the inter-dot coupling 
is well approximated by the spin coupling of two dots. 

In Appendix A, we demonstrate the absence of the 
pumping current in the single quantum dot, and the va- 
lidity of the Kondo model when the temperature is much 
higher than the characteristic energy, so called Kondo 
temperature. 

Our model is depicted in Fig. [TJ The left and right 
quantum dots are connected with the lead L and lead R 
by the spin couplings Jl and Jr, respectively. Here we 
ignore the energy dependence of Jl and Jr by renormal- 
izingthe high energy fluctuation into the couplings (see 
Ref . l39i ) . Thus Hamiltonian is given by 

H = H a + Hkl + Hkl + H12, (1) 

H = e k a \,k,a a l,k,<7, (2) 

-Hkq = Jl&1 ■ SL,k,k' + JrS2 • Sfi^fc/, (3) 

k,k' 

H12 = JiaSx-Sa, (4) 

where a) k a ,a 7 ^,a are the creation and annihilation op- 
erator of the electron with the wave number k and the 
spin a in the lead j{= L,R), Si (S2) is the spin opera- 
tor in the dot 1 (2), and sx,k,k' = J2 S , S < a l.k,s a s,s'a 7 ,k,s' ■ 
Here we adopt the Hund coupling as the interaction be- 
tween the quantum dots. The derivation of the Hamilto- 
nian (jl])- (J^J) is presented in Appendix. 

We consider an adiabatic modulation of the param- 
eters in the outer region of the double quantum dots. 
Throughout this paper, we treat the temperature in each 
leads (Tl and Tr) and chemical potentials (/ij, and /j,r) 
as the control parameters. It is straightforward to expand 
the parameter space to include the parameters in inner 
region of the quantum dots, e.g., energy levels, intra-dot 
couplings, tunnel potential that corresponds to the mod- 
ulation of the spin couplings. 



FIG. 1: A schematic picture of double quantum dots coupled 
to external leads 
A double quantum dots with spin 1/2 in each is connected 
to the external leads. The inter and intra dot interactions 
are considered as the spin coupling. We have two control 
parameters in each lead; temperature and chemical 
potential. 

III. METHOD 

In this section, we make a brief review of the method 
to analyze the geometrical pumping in the context of 
QME.— Here we apply the Keldysh method for calcula- 
tion of the generating functionj 33 i 34 We treat the left and 
right leads as the heat and particle baths. After repre- 
senting the generating function by the path integral, we 
integrate out the Grassmann variables for the leads. The 
general results in the following coincides to that used in 
Refs. [H and [3^. However, that is more tractable and 
systematic than the method used in Ref. [H. 

First, we count the number of the electrons transfer 
from the left lead to the right lead by introducing the 
counting field. The generation function is given by 

Z(x) = ( e - lxA V xJVo ), (5) 

where Nq and N t are the number of the electrons in left 
leads at initial time and at time i, respectively. For our 
model, the commutable relation [Ho,N] = holds. Thus 
Eq. ([5]) is rewritten as 

Z( X ) = Tr [e^e-^Vo] , (6) 
where po is the density matrix for initial state and 

e iH ±x t _ e ±ixN e iHt e T± x N _ ( 7 ) 

Equation ([5]) denotes the time evolution depicted in Fig. 
[2] The initial state described by p evolves by the 
"Hamiltonian" H_ x from initial time to t and the state 
at t evolves backward by H x . By using the Keldysh's 
method of contour integral, we obtain 

Z( X ) = Tv s T c J DCD^p 0s e-^ (G ' 1+H ^, (8) 

where the contour c is represented in Fig. [2l T c is the 
contour ordering operator defined on c. The bare Green's 
function is defined by 

G = <a'(i9t-H )-\ (9) 
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IV. GEOMETRICAL PUMPING 
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In this section, we apply the method described in the 
time previous section to the model depicted in Fig. 1. In the 

first subsection, we derive the QME for our system. We 
H-^ apply the Markovian approximation and perturbation for 

spin coupling, we obtain the QME with counting field. In 
the second subsection, we calculate the BSN phase and 
FIG. 2: The contour of integral appears in the generating curvature for our system. We also estimate the pumping 
function (0. current for cyclic modulation. 



where a z represent the z component of the Pauri matrix 
and a, a' denotes the upper contour or lower contour. By 
integrating out the Grassmann field, we obtain 



Z( X ) = Tr s T c p 0s e Trln ^ 1+ ^. 



(10) 



Here we can define the density matrix for the system as 
p x = T , c pose Trln ( G +h k). By using the Markov approx- 
imation, the time evolution of p x becomes local in time 
and thus we obtain the following formal quantum master 
equation 



dt s 



(11) 



Next we consider the modulation of the parameters. 
We follow the method described in Ref. The density 
matrix for the system can be expanded as 



P x (t) 



(12) 



where at denotes the set of parameters at time t, X n and 
p x is the eigenvalue and corresponding eigenfunction of 
the K, in Eq. $TB). Here, A*(t) = f*dt'\ x (a t ) is the 
dynamical phase corresponds to the house-keeping part 
of the cumelant generating function. Thus we obtain the 
excess part by subtracting the house-keeping part from 
the total cumelant generating function. The excess part 
of the cumelant generating function becomes 

Sex(ix) = - / \ l o,s(t)' d Po,S(t)) + surface terms, (13) 



where p x s ^ and l x s , t s are the right and left eigenvectors 
corresponding to the maximum eigenvalue, respectively. 
Here a(t) is the vector which consist of the control pa- 
rameters. Thus we obtain the average number of elec- 
trons transfered from the left lead to the system AN by 
differentiating Eq. ([TB")) with respect to i\, 



AN : 



(^0,<3(t)' ^Po,S(t] 



(14) 



where l 'o,s(t) = dt o,sy)/ d ( i x)\x=o- 



A. Quantum Master Equation 

We introduce the counting fields x 7 o- to count the num- 
ber of the electrons N ja - with spin a in lead 7. The gen- 
erating function becomes 

Z( X ) = K s T cPOs e^ G - 1+H l« +H l« +H ^, 
with 



(15) 



R x - T • S S ^ a n z 



i(Xio 



,)o 



>/2 



(16) 



where a = + or — represents the Keldysh branch and 

S^J?" = Jdte-^-^Si^it). Here we define 5{L) = 
1, 8{R) = 2. The Hund coupling term H±2 is commutable 
with the number operator for leads. We focus on the case 
of ttvJj <C 1 and JlrJ J-y <C 1. Thus we make the per- 



turbation of J 7 and J12 for Eq. (|15|) . We can show that 
the off diagonal part of the density matrix is irrelevant in 
the steady state. In the case of adiabatic modulation we 
have assumed, the contribution of these terms is negligi- 
ble. Thus the QME is closed in the diagonal part of the 
density matrix, p x a , , where a and a' stand for the spin in 
the left and right dots, respectively. By using the Marko- 
vian approximation, the QME for the diagonal part of p x 
becomes 



d t p x s = R x p x 



R x 



(a 




bt 





\ 


b~ R 


a 


—Jlr 


bl 






—Jlr 


a 


b + R 




\ 


K 


bR 


a 


) 



(17) 



(18) 



in the leading order of the perturbation, where, p x — 

(ptt>Pn>Ptt>Pu)* and 



7,7' 

x jYg-^Xrl— X 7 'f) -|- 6 -*(Xt4.-X7'4.)^ _ g 



, (19) 



de 



= (^J 7 ) 2 J — [1 - / 7t (e)] f^(e)e-^- x ^0) 
= (™J 7 ) 2 /"^[l-/ 7 4(e)]/ 7t (e)e-^-H ? i) 
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with Fermi 

[l — e -/3 7 (e-A*T)" 



distribution 

-l 



function / 7 (e) 



T -i 

7 



in the lead 7. Here j3 1 
and fi~[ are the inverse temperature and the chemical 
potential in the lead 7, respectively. As we men- 
tioned above, here we consider the parameters to be 
a(t) = (pl,T l ,h r ,T r ). 



B. BSN curvature and pump current 




It is easy to show that the QME (jl8|) has four eigen- 
values A a (a = 0, 1, 2, 3), 



Ao = 
Ai = 



2a- 


K&Z- 


!-#) + 




-b + L f- 


-4^6+ 


,(22) 


2a - 


"(61- 






-b + L y- 


-±b- R b + R 


,(23) 


2a- 


Hbl- 


vb + L )- 




-b + L y- 


-±b~ R b + R 


,(24) 


2a 


-<PI 






-bir- 


^b R b+ 


.(25) 



A 2 = - 



A, = - 



It is obvious that the eigenvalue Ao has the largest real 
part. When X = 0, the right eigenstates corresponding 
to Ao becomes, 

Po = -(1 + A,1-A,1-A,1 + A)\ (26) 

where, A = 2 Jlr/(1>l + b R ) is dimensionless parame- 
ter with 6l = bf^\ x= o,b R = o^| x= o and t represents the 
transverse. The quantities 6^ and b R can be calculated 
by Eq. {H as 



T 

& 7 = (vri/J 7 ) 2 ^, 



(27) 



The derivative of left eigenvalue with respect to ixLf at 
X = which enters in Eq. (fbfj) becomes, 



ft = (-^,o,Ko 



fe-fe-v (28) 

The BSN curvature for our model formally satisfies 

5^ i^dz A dy = dl' A dp , (29) 

where x,y — /il, Pr, Pl> Pr, and the total derivative in 
the right hand side is taken in the parameter space and 
A is the edge product. Substituting Eqs. (f26|) - (f28)) into 
Eq. (|2"9"]l . the BSN curvature becomes 



Tl+P 2 lr T r 



(30) 



where a±2 — J12/ (t^^Jr) 2 and (3lr = Jr/Jl- This ex- 
pression of the curvature is proportional to Tj, at high 
Ti,. It is notable that the curvature always survives 



FIG. 3: The density plot of curvature Fp L p R /a\2 in Eq. (f30|) 
in a case of Jl = Jr- The curvature is nonzero in the tem- 
perature region where the temperature is much higher than 
the Kondo temperature. 



thanks to the spin coupling. We should note, however, 
that our model is valid for spin localized region i.e. for 
Tl,T r <C U, where U is the intra-dot interaction. Thus 
we expect that the Eq. (|3"0|) saturates to finite value. Our 
model in which the two Kondo impurities are coupled 
with spin coupling requires the condition J 12 <C Jl, Jr. 
In Fig.|3l we show the curvature Eq. (|30[) . The curvature 
is non-zero in the whole parameter region where the tem- 
perature is much higher than the Kondo temperature. 

In our system, the BSN curvature vanishes in the chem- 
ical potential plane. This result seems completely differ- 
ent from the results in Ref. [32T which offers non-zero 
curvature in chemical potential plane. However, the re- 
sult we obtain above is consistent with the previous re- 
sults by following reason. The BSN curvature in Ref. [32| 
has peaks around energy levels 61,62 in quantum dots 
and ei + U, 62 + U. The result clearly shows that the 
pumping is assisted by the charge fluctuation. In this 
paper, however, we restrict our model to the spin lo- 
calized region by suppressing the charge fluctuation in 
construction of Kondo model. Thus our model does not 
allow us to describe the quantum pumping assisted by 
charge fluctuation which is dominant in chemical poten- 
tial modulation. In our case, the spin fluctuation assists 
the quantum pumping. 



C. Pumping current 

In this subsection, we explicitly calculate the pumped 
electrons in the left lead, ANl^ generated from t = 
to to. In order to calculate AN L ^ we set xl-[ = X and 
Xli = 0, XRcr = as in the previous section. Thus 
we obtain l' Q = (— b R /bL, 0, b R jb^, 0). Substituting the 
above expression of 1' into Eqs. ([2"o) - (|2"7) into Eq. fiT 
we obtain 



AiV it = J L r J ^d(b L 



b R y 



(31) 



As we show in the previous subsection, the BSN cur- 
vature in the parameter space is non-zero. Thus ANl^ 



AA/f 




FIG. 4: The excess current for square modulation in temper- 
ature plane of the parameter space. The left figure shows 
the trajectory in the temperature plane. We start from 
(T L ,T R ) = {T , T ) and increase T R from T to T x . Next we 
change Tl as To —¥ T\. Then we modulate the temperature 
from (Ti, Ti) to (To, Ti) by lowering Tl and go back to origi- 
nal point by lowering Tr. Here we set Jl = Jr, To = 0.5 and 
ai2/2To = 1. The right figure shows the Ti dependence of 
the excess current. We show that the larger area enclosed by 
the trajectory achieves larger creation of the excess current. 

obtained by ([3"Tjl depends on the trajectory in the pa- 
rameter space. As an example, we calculate the excess 
current for modulation depicted in left figure of Fig. [4] 
We change the temperature in left and right leads from 
(T , T ) to (Ti, Ti) by the square trajectory in Fig.|U In 
the right figure of Fig. 0] we show the Ti dependence of 
the excess current achieved by the one cycle. The re- 
sults shows that we can obtain the large excess current 
by raising T\. Here we again mention that the range of 
the temperatures To and Ti are restricted in the region 
much larger than Tk and much less than U . 



V. CONCLUSION 

We have derived the geometrical pumping for a sys- 
tem of double quantum dots connected with heat baths 
by external leads based on the QME approach. In the 
spin localized region, the spin fluctuation is dominant for 
the quantum pumping and that is complementary to the 
result for a spin less model in Ref. [H where the charge 
fluctuation is dominant. We have also obtained the an- 
alytical expression of BSN curvature in the parameter 
space for our model. We have shown that the geometrical 
pumping is not realized by the single dot nor the isolated 
two dots (double dot without inter dot repulsion) by the 
modulation of temperature or chemical potential, where 
each dot is restricted to the spin localized region and has 
spin degeneracy. According to our analysis, if we destroy 
the spin degeneracy by applying the magnetic field or 
connected to the ferromagnetic leads, the spin pump can 
be yielded. In contrast to that, we obtain the pumping 
induced by the inter dot repulsion. In this case, we ob- 
tain non zero pumping by modulation of the temperature 
in the external leads. 

In this paper we restrict our analysis to the quan- 



5 



dot 




Jl Jr 

FIG. 5: Model for the single dot connected to the leads. 

turn dots in the spin localized region (weak coupling re- 
gion). However, the QME formalism is applicable to the 
whole region of the coupling strength. In the strong cou- 
pling region, the Kondo singlet state is formed around 
the Fermi surface, that is not captured by the spinless 
model in Ref. [32| nor the weak coupling model in our pa- 
per. The Kondo singlet state is well known to bring a 
drastic change to transport phenomena at the tempera- 
tures lower than the Kondo temperature. Concerning to 
that, the spin pumping in the Kondo regime for Toulouse 
limit has been proposed recently. 17 It is a future issue to 
reveal the effect of the many body correlation on the 
quantum pump. The investigation of the non-adiabatic 
correction and non-Markovian effect are also the future 
issues. In the case of boson transport, the non-adiabatic 
effect and non-Markovian effect have been studied and 
the generalization to the case of fermion transport is 
straightforward^ 
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Appendix A: Case of single dot 

We devote this appendix to show that the geometrical 
pumping is prohibited for the single quantum dot system. 
Our model analyzed in this appendix is depicted in Fig. 

m 

The Hamiltonian is given by 
H = H a + H int , (Al) 

H T = ^2v a 4a L + h.c. , (A3) 

a, cr 

where, a} aka and a a .k, a are the creation and the anni- 
hilation operator for the electron in the leads a(= L,R) 
with wave number k, energy ej,, and spin a. d\,d a is 
those in the quantum dot and n a — d\d a . U is the 
electron-electron interaction in the quantum dot. 
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By using the Schriefer- Wolff transformation 3 ^, 
above model reduces to the following model: 

H = H n + Hk, 

Hq = ^ £k a k,a a k,cn 
k.a 



Hk — ^ JS ■ Sfc,, 



the 

(A4) 
(A5) 

(A6) 



k,k' 



where, a kA 



V L a L „ 



Vror^/V, V 



\J\Vl\ 2 + \Vr\ 2 and S is the spin operator in the quan- 
tum dot. We denote the spin operator for the conduction 

electrons as s k:k , = J2 a ,a> <*k,<j a W ,°' a ",v' l^- Tne coupling 
constant J is presented by the density of states in leads v 
as J = lirvV 2 1 (eo + U). We assume the density of states 
to be constant (wide-band limit). Here we focus on the 
electron- hole symmetry (|eo| = eo + U) for simplicity. 37 

First, we renormalize the spin fluctuation from the 
high energy region. We apply the poor man's scaling 
method 3 ^^ to the Kondo model Eq. (|A"4)) . This yields 
the effective coupling 



J 



2i/ mmax{e V,T}/T K ' 



where V is the bias voltage between the left lead and the 
right lead. The Kondo temperature Tk is defined by 



T K = D exp 



1 

2^7 



(A8) 



In this paper, we focus on the spin localized region (or 
the weak coupling region) vJ < 1. In the region, the 
conduction electrons are scattered by the local spin in 
the dot sequentially. 

In this case, the QME for the diagonal parts of p x 
becomes 



a (■KvJ) 2 b+ 
{irvJ) 2 b~ a 



(A10) 



in the leading order of the perturbation in J by using the 
Markov approximation. Here a is the given by Eq. f| 19[) 
with Jl — Jr — J and 







7,7 




7,7 





b~ = 



Thus we obtain the right and left cigcnstates correspond- 
ing to the eigenvalue with largest real part as 



Ps(t) 



lf^/b+Jb- 



%) = ( Vb=W 1 ) 



(A13) 
(A14) 



(A7) Taking the derivative with parameters for p%i t \ yields 



d Ps(t) 



1 



d3.| m(VbV^) 

V2 \ 



(A15) 



By setting \ — and using the property 6 + | x= o 
b~\ x= o, we obtain 



dpi 



S(t) 



= 0. 



(A16) 



dtp* = R*ri 



(A9) Thus the curvature vanishes in the parameter space. 
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